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Abstract 

New ^^CO J=4-3 and ^^CO J=3-2 observations of the N159 region, an active 
site of massive star formation in the Large Magellanic Cloud, have been made with 
the NANTEN2 and ASTE sub- mm telescopes, respectively. The ^^CO J =4-3 dis- 
tribution is separated into three clumps, each associated with N159W, N159E and 
N159S. These new measurements toward the three clumps are used in coupled calcu- 
lations of molecular rotational excitation and line radiation transfer, along with other 
transitions of the ^^CO J=l-0, J=2-l, J=3-2, and J=7-6 as well as the isotope 
transitions of ^^CO J=l-0, J=2-l, J=3-2, and J=4-3. The ^^CO J=3-2 data are 
newly taken for the present work. The temperatures and densities are determined to 
be ~ 70-80K and ~ 3 x 10^ cm'^ in N159W and N159E and ~ 30K and ~ 1.6 x 10^ 
cm-3 in N159S. Observed ^^CO J=2-l and ^^CO J=l-0 intensities toward N159W 
and N159E are weaker than expected from calculations of uniform temperature and 
density, suggesting that low-excitation foreground gas causes self-absorption. These 
results are compared with the star formation activity based on the data of young 
stellar clusters and H II regions as well as the mid-infrared emission obtained with the 
Spitzer MIPS. The N159E clump is associated with embedded cluster(s) as observed 
at 24 fim by the Spitzer MIPS and the derived high temperature, 80K, is explained 
as due to the heating by these sources. The N159E clump is likely responsible for 
a dark lane in a large Hll region by the dust extinction. On the other hand, the 
N159W clump is associated with embedded clusters mainly toward the eastern edge 
of the clump only. These clusters show offsets of 20" - 40" from the ^^CO J=4-3 peak 
and are probably responsible for heating indicated by the derived high temperature, 
70 K. The N159W clump exhibits no sign of star formation toward the ^^CO J=4-3 
peak position and its western region that shows enhanced i?4_3/i_o and -R3_2/i-o ra- 
tios. We therefore suggest that the N159W peak represents a pre-star-cluster core of 
~ IO^Mq which deserves further detailed studies. We note that recent star formation 
took place between N159W and N159E as indicated by several star clusters and Hll 
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regions, while the natal molecular gas toward the stars have already been dissipated 
by the ionization and stellar winds of the OB stars. The N159S clump shows little 
sign of star formation as is consistent with the lower temperature, 30K, and some- 
what lower density than N159W and N159E. The N159S clump is also a candidate 
for future star formation. 

Key words: Magellanic Clouds — ISM: clouds — ISM: molecules — radio lines: 
ISM — submillimeter — stars: formation 

1. Introduction 

Giant molecular clouds (hereafter GMCs) are the principal sites of star formation and 
studies of GMCs are important in understanding the evolution of galaxies. A few tens of GMCs 
in the solar neighborhood such as Orion A have been well studied, while most of the GMCs in 
the Galaxy are located at more than a few kpc away in the Galactic disk, where contamination 
by un-related components in the same line of sight seriously limits detailed understanding of 
GMCs and their associated objects. 

Recent CO surveys of molecular clouds toward external galaxies in the Local Group have 
revealed that properties of GMCs such as relations between CO luminosity and line width, the 
mass range, and the index of mass spectrum are similar among these galaxies (Blitz et al. 
2007). This suggests that the properties of GMCs are fairly common among galaxies. The 
Large Magellanic Cloud (hereafter LMC) is the nearest galaxy to our own at a distance of 50 
kpc and is nearly face on, making it an ideal laboratory to observe various properties of GMCs. 
The small distance enables us to make molecular observations at high spatial resolutions. The 
LMC also provides a unique opportunity to study molecular clouds and star formation in 
different environments from the Galaxy; the gas-to-dust ratio is ~ 4 times higher (Koornneef 
1982), and the metal abundance is about ~3 - 4 times lower (Rolleston et al. 2002; Dufour 
1984) than those of the Galaxy. 

Following a low resolution ^^CO J=l-0 survey at 140 pc resolution with the Columbia 
1.2 m telescope located at CTIO in Chile (Cohen et al. 1988), high resolution observations in 
the ^^CO J=l-0, J=2-l emission line were made with the SEST 15 m telescope toward some 
of the molecular clouds and revealed their clumpy structure at 10 pc resolution (e.g., Israel 
et al. 1986; Johansson et al. 1994; Caldwell & Kutner 1996; Kutner et al. 1997; Johansson et al. 
1998; Israel et al. 2003). However, these high resolution studies were limited in the spatial 
coverage, compared with the large angular extent of the LMC, ~ 6° x 6°. Fukui et al. (2008) 
carried out a survey in the ^^CO J=l-0 emission line at 40 pc resolution over a 6° x 6° field in 
the LMC with the NANTEN 4 m telescope and obtained a complete sample of 270 GMCs (see 
also for preceding works Fukui et al. 1999; Fukui et al. 2001; Mizuno et al. 2001; Yamaguchi 
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et al. 2001). These studies revealed that young stellar clusters whose ages are less than 10 Myr 
are spatially well correlated with GMCs, and that GMCs are categorized into three types in 
terms of star formation activity. Type I is starless in the sense that they are not associated 
with O stars, Type II is associated with small H II regions only, and Type III is associated with 
huge Hll regions and young stellar clusters (Kawamura et al. 2009; Blitz et al. 2007). These 
types are interpreted in terms of evolutionary sequence of GMCs in a timescale of 2-30 Myrs 
(Fukui et al. 1999; Blitz et al. 2007; Kawamura et al. 2009). 

The ^^CO J=l-0 emission line is a probe commonly used to trace molecular clouds 
because of its low excitation energy (~ 5K) and low critical density for collisional excitation 
(ncr ~ lOOOcm"'^). The ^^CO J=l-0 emission alone is however not able to provide physical 
properties like kinetic temperature and density, the fundamental parameters of GMCs. The 
high J transitions have higher excitation energies and higher critical densities; e.g., the ^^CO 
J=3-2 transition has the upper state at 33 K and the critical density, 3 x lO'^cm"^ (Schoier 
et al. 2005), and the ^^CO J=4-3 transition has the upper state at 55 K and the critical 
density, 1 x lO^cm"'^. These sub-mm CO emission lines can selectively trace the sites which 
may be warmer and denser than the mm CO lines and have a potential to reveal detailed 
physical properties where star formation is taking place. We are also allowed to make better 
estimates of temperatures and densities of molecular clouds with a combination of multi-J CO 
line intensities and molecular excitation analyses such as the Large Velocity Gradient (hereafter 
LVG) model of molecular line transfer. Recently-developed sub-mm telescopes such as ASTE, 
NANTEN2 and APEX located at altitudes around 5000m in Atacama, Chile have enabled 
us to observe higher excited CO transitions at sub-mm wavelengths in superb observational 
conditions (Ezawa et al. 2004; Fukui et al. 2006; Giisten et al. 2006). At a lower angular 
resolution, AST/RO also offered sub-mm observing capability (Stark et al. 2001). 

Bolatto et al. (2000) first detected the ^^CO J=4-3 emission line toward the N159 region 
at 50 pc resolution with the AST/RO telescope. Bolatto et al. (2005) later presented estimates 
of temperatures and densities; T^in = 20 K, n(H2) = 10^ cm^'^ for the cold dense component 
and Tkin = 100 K, n(H2) = 100 cm~^ for hot tenuous component toward N159W. Subsequently, 
Minamidani et al. (2008) carried out high resolution ^^CO J=3-2 observations of several GMCs 
at 5 pc resolution with the ASTE 10 m telescope including N159 in the LMC. Minamidani et al. 
(2008) revealed detailed structure of highly excited gas with Tkm ~ 20 - 200 K and n(Il2) ~ 
10^ - 10^ cm~^. In the present study, we shall focus on the N159 CMC which shows the 
highest ^^CO J=l-0 intensity from the NANTEN survey (Fukui et al. 2008). The N159 CMC 
is classified as Type III, and includes at least three prominent clumps, N159W, N159E and 
N159S (Johansson et al. 1994). The two molecular clumps in the northern part, N159W and 
N159E, are associated with massive star clusters whose ages are younger than 10 Myr (Bica 
et al. 1996) and with huge Hll regions (Henize 1956; Davies et al. 1976; Kennicutt & Hodge 
1986). On the other hand, there is no star formation in N159S. In order to better estimate 
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the physical properties of the N159 region, we have carried out new sub-mm observations in 
the ^^CO J=4-3 emission hne at a 10 pc spatial resolution with NANTEN2 and in the ^^CO 
J=3-2 emission line at a 5 pc spatial resolution with ASTE. We have also made combined 
calculations of molecular rotational excitation and line radiative transfer by employing the CO 
datasets available in N159 and compared the results with the star formation activity. 

The present paper is organized as follows: Section 2 describes the observations. Sections 
3 and 4 show the observational results and data analysis, respectively. We discuss the correlation 
between highly excited molecular gas and star formation activities in section 5. Finally, we 
present a summary in section 6. 

2. Observations 

2.1. ^"^ CO J =4-3 observations 

We observed the ^^CO J=4-3 rotational line at 461.0408 GHz with the NANTEN2 
(Fukui et al. 2006) 4 m sub-mm telescope situated at a 4800 m altitude at Pampa la Bola in 
Chile. The telescope has a main beam efficiency of 0.6 and the half power beam width (HPBW) 
of 38", as determined from continuum cross scans on Jupiter at 461 GHz. The forward beam 
coupling efficiency was measured by sky-dip measurements to be 0.86. We used a dual-channel 
460 / 810 GHz receiver developed by University of Cologne. Observations were made during 
September to December in 2006 for N159W and October to December in 2007 for the N159E 
and N159S regions. The double sideband receiver temperature was measured to be 250 K at 
460 GHz and the typical system noise temperature during these observations was ~ 900 K 
toward the zenith. The spectrometer was an acousto-optical spectrometer (AOS) with total 
bandwidth of 1 GHz divided into 1730 channels, yielding a velocity resolution of 0.37 km s~^. 
We observed three 2'x2' areas toward the three prominent clumps; N159W, N159E and N159S 
(Johansson et al. 1998). Observations were made using on-the-fly (OTP) mode with two 
orthogonal scan directions along RA and DEC to ensure a flat fleld and the data were gridded 
with a 10" spacing. Calibration and data gridding were made with the "kalibrate" software 
developed by University of Cologne. The pointing accuracy was better than 10", as conflrmed 
from pointing measurements on IRC2 in Orion A. The flnal noise levels are 1.2 K/ch for N159W 
and 0.4 K/ch for N159E and N159S. 

2.2. ^^CO J =3-2 Observations 

The ^'^CO J=3-2 rotational line at 330.5879 GHz was observed with the ASTE 10 m sub- 
mm telescope in Pampa la Bola in Chile (Ezawa et al. 2004; Ezawa et al. 2008), using a double 
side band superconductor insulator superconductor (SIS) mixer receiver (Muraoka et al. 2007) 
in a period from 20 to 23 September 2006. The half-power beam width was measured from 
observations of Jupiter to be 23". The spectrometer was an XF-type digital auto correlator 
(Sorai et al. 2000) and was operated in the wideband mode, having a bandwidth of 512 MHz 
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with 1024 channels. This configuration yielded a velocity coverage and resolution of 450 and 
0.45 km s~^, respectively. We observed 3x3 positions at 20" spacing in the position switching 
mode centered at the ^^CO J=3-2 integrated intensity peaks (Minamidani et al. 2008) in 
N159W (5^39'°36.8% -69° 45' 32" at J2000 coordinate), N159E (5'^40°'8.7^ -69° 44' 34") and 
N159S (5M0™5^, -69° 50' 34"), where the observing grid is along the B1950 coordinate. The 
pointing error was measured to be within 5" from the observations of the ^^CO J=3-2 point 
source R-Dor. Typical system noise temperature was measured to be 260 K in DSB. Finally, the 
map was re-gridded to the J2000 coordinate by using "miriad" astronomical data analyzing 
software. We observed Ori-KL and M17SW to check the stability of the intensity calibration, 
and the maximum intensity variation during these observations was 12%. We use 0.7 for main 
beam efficiency, as measured from observations of Jupiter. Due to a lack of side-band ratio 
data, the calibrated intensity scale was somewhat uncertain in 2006, and we re-observed the 
peak position of N159W with the newly implemented single side band receiver (Ezawa et al. 
2008) in 2008, and estimated the main beam intensity to be 3.15 ± 0.15 (K). We scaled the 
observational data in 2006 by using the derived scaling factor. Achieved noise levels in the 
main beam temperature scale are 0.5 K/ch for N159W and 0.2 K/ch for N159E and N159S. 

3. Observational Results 

3.1. The ^^CO J =4-3 and ^^CO J =3-2 emission lines 

Figure 1 shows the spatial distributions of the ^^CO emission lines. The right panel 
of Figure 1 shows the distribution of the ^^CO J=4-3 integrated intensity smoothed to a 40" 
beam, assuming a Gaussian beam shape. For comparison, the left and middle panels of Figure 1 
show the distributions of the total integrated intensity of the ^^CO J=l-0 (Ott et al. 2008) and 
the ^^CO J=3-2 (Minamidani et al. 2008), respectively. Distributions of these ^^CO emission 
are largely similar with each other. The ^^CO J=4-3 distribution has clearly resolved the 
cloud into three clumps, N159W, N159E, and N159S which were identified in the ^^CO J=l- 
observations (Johansson et al. 1994). These three were not resolved in the previous ^^CO 
J=4-3 observations (Bolatto et al. 2000). 

The ^^CO J=3-2 emission was measured in a 3 x 3 pointing pattern, with a 20" grid 
spacing toward each ^^CO J=3-2 peak of N159W, N159E and N159S. The line profiles are 
shown in Figure 2. 

Figure 3 shows four spectra of the ^^CO J=l-0, J=3-2 and J=4-3 and ^^CO J=3-2 
emission lines toward three peaks of the ^^CO J=3-2 integrated intensity (Minamidani et al. 
2008). The spectra are convolved to a 45" Gaussian beam for a quantitative comparison. 
N159W shows the strongest emission among the four spectra. 

Table 1 lists the observed parameters of the ^^CO J=4-3 and ^^CO J=3-2 emission 
lines toward selected peak positions of the N159 region. The parameters of the other emission 
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lines including ^^CO J=l-0, J=2-l, J=3-2 and J=7-6 and ^^CO J=l-0, J=2-l, and J=4-3 
(Johansson et al. 1998; Pineda et al. 2008; Ott et al. 2008; Minamidani et al. 2008) are also 
summarized in Table 1. Nine transitions are available in N159W and five in N159E and N159S. 

4. Data Analysis 

4.I. Large Velocity Gradient Analysis 

In order to estimate physical properties of molecular gas, we employ a LVG model 
(Goldreich & Kwan 1974). A LVG model is useful to calculate level populations of quantum 
states of a molecule and line intensities of each transition. A model molecular cloud is assumed 
to have a large velocity gradient as compared with local velocity dispersion. The line emission 
from one place affects via absorption or induced emission only the local region whose size is 
in the order of the local velocity dispersion divided by the velocity gradient. By using an 
escape probability which is the probability for a once-emitted photon to escape from the cloud 
without re-absorption, the equations of the line radiative transfer and the statistical equilibrium 
are simplified and are solved locally. It is then possible to calculate the population of excitation 
states. 

Here, we calculate level populations of ^^CO and ^^CO molecular rotational quantum 
states and line intensities under an assumption of a spherically symmetric molecular cloud of 
uniform density and temperature. The escape probability for a spherically symmetric velocity 
gradient is calculated as 

1 - exp{-Tj+ij) 

Pj+i,j = i-Lj 

TJ+1,J 

(Castor 1970). The LVG model requires 3 independent physical parameters to calculate emis- 
sion line intensities; kinetic temperature, gas density and X/{dv/dr). X/{dv/dr) is the molec- 
ular abundance divided by the velocity gradient in the cloud; [^"^ CO / B.2] / (dv/dr) for ^^CO and 
(p2CO/H2]/p2co/i3CO])/(di;/dr) for ^^CO. 

We cannot measure molecular abundances directly. Here, we estimate [^^C0/H2] from 
the Aco-factor derived by Fukui et al. (2008). The Aco-factor in the LMC is about 3 times 
higher than the Galactic value and we shall adopt [^^C0/H2] in the LMC 3 times lower than the 
Galactic value. Since the Galactic Aco value was determined in the Orion molecular cloud to 
be 5 X 10"^ by Blake et al. (1987), the estimate [^^cO/Hs] in the LMC is 1.6 x 10"^ Johansson 
et al. (1994) estimated that the [^^CO/^^CO] value in the N159 region is about 2/3 of the 
Galactic value from observations of the rare CO isotopes (C^^O, C^^O). Thus, by assuming 70 
for the Galactic value we adopt [^^CO/^^CO] to be 50 in the N159 region. The velocity gradient 
within the molecular clumps in the N159 region is typically 0.5 km s"^ pc^^ (Minamidani et al. 
2008). Therefore, we adopt X/{dv/dr) = 3.2 x 10"^ for ^^CO and 6.4 x 10^^ for ^^CO. 

In order to obtain the optimum physical parameters, we calculate values from observed 
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intensity ratios derived from intensities which are convolved to the same beam size (45") at 
peak positions of the ^^CO J=3-2 emission (Table 2). 



where N is the number of observed molecular transitions and "i" and "j" refer to different molec- 
ular transitions. Robs{i,j) is the observed line intensity ratio of transition "i" over transition 
"j", and Rivdhj) is the ratio of the LVG calculations at any temperature and density, "a" 
is the uncertainty of the observed line intensity ratio. The errors in the observed intensity are 
mainly due to calibration uncertainty; we use conservatively a as 20% of the peak line intensity 
ratio. Table 1 shows observed line intensity at a peak position of each clump that was used in 
the LVG calculations. All intensities were smoothed to a 45" beam size. 

4-2. Results of the LVG Analysis 

We calculated the reduced using all observed line intensities listed in Table 1. We 
define the range of solutions to be within the 1 a level of the reduced value. We display the 
reduced distributions in gray scale and the range of solutions in red contours in the density- 
temperature plane in Figure 4 (left). A subset of observed line intensity ratios is summarized 
in Table 2 and also plotted in Figure 4 (green lines). From the minimum of reduced values 
in each panel, we estimate density and temperature at the three peaks as given in Table 3 and 
find that N159W and N159E {Tkin ~ 70 - 80 K and n(H2) ~ 3 - 4 xlO^ cm'^) are significantly 
warmer than N159S {Tkm ~ 30 K and n(H2) ~ 1.6 xlO"^ cm~^). The reduced values of ~1 
for N159W and N159E indicate that the solutions are not very reliable. 

We note that locus C, the ratio between the ^^CO J=4-3 and J=3-2 lines, deviates from 
the solution in M159W (Figures 4 a and b). We suggest that this can be a result of calibration 
error of 20% (see Table 2). We also note that the inclusion of loci A and B in Figures 4 a, c, 
and e, the ratios including the ^^CO J=l-0 and J=2-l transitions, tend to make the fitting 
somewhat worse (i.e., = 1 in N159W and N159E). We repeated the calculations by excluding 
intensity ratios with ^^CO J=l-0 and J=2-l. The solutions are shown in Figures 4 b, d, and f. 
Thus obtained temperatures and densities are summarized in Table 3. The results are similar to 
those obtained for the original calculations, but the reduced x^ value is improved (i.e. x^ = 0.5 
in N159W and = 0.07 in N159E). We suggest that the lower excitation mm emission lines 
may be affected by the foreground gas different from that emitting the sub-mm emission lines 
as shown later in Discussion. 

We shall then discuss on the possible effects in determining temperature and density 
due to the difference in observed positions among different transitions (Table 2). The observed 
positions in ^^CO J=2-l, ^^CO J=l-0, and ^^CO J=2-l by Johansson et al. (1998) are shifted 
from those in the present work by ~10 " in N159W and N159E and by ~30 " in N159S for 




{RpbsjiJ) - RLVG{iJ)V 

a2 



(2) 



i=l j=i+l 



Reduced x^ = X^/{N-l) 



(3) 
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the 45 " convolved beam. The intensity differences each between these two positions are less 
than 5% according to the OTF distributions of the ^^CO J=l-0, J=3-2, and J=4-3 lines. In 
LVG calculations we find the resultant differences in temperature to be ~5% in N159W and to 
be ~10% in N159E and N159S and that density remains the same, if the ^^CO J=2-l, ^^CO 
J=l-0 and ^'^CO J=2-l line intensities are shifted by 10% from the values of Johansson et al. 
(1998). We summarize that the effects of the positional differences in Table 2 may cause errors 
up to 10% in temperature, while no significant effect in density. 

4-3. Discussions on the LVG Results 

We show a comparison of the calculated and observed intensities as a function of J in 
Figure 5. This comparison indicates that the line intensities in the high-excitation sub-mm 
regime are well reproduced and we argue that a combination of ^^CO J=3-2, J=4-3, and ^^CO 
emission lines are crucial to derive physical properties of the embedded dense region. On the 
other hand, in N159W and N159E the lower transitions, ^^CO J=l-0 and J=2-l, are measured 
to be significantly below the calculated values. The ^^CO J=l-0 line of N159W in Figure 
3 shows a somewhat fiat-top profile, suggesting possible self absorption by some foreground 
component which may cause the lower intensity of the ^^CO J=l-0 emission in N159W and 
N159E (Figure 5), while we admit the noise levels are higher in the ^^CO J=l-0 spectra than 
the others. For instance, LVG calculations for density of 100 cm~^ and temperature of 30K 
with the same set of the other model parameters indicate that the 20 % decrease of the J=l-0 
and J=2-l lines in N159W is well reproduced by self-absorption by such diffuse gas, while the 
density and temperature of the foreground gas is not necessarily uniquely determined. Such 
foreground gas is not significant in emitting the J=3-2 and J=4-3 lines. 

In N159S, on the other hand, all the observed intensities are consistent with calculations. 

We thus conclude that the N159W and N159E clumps consist of warm and dense gas 
with a temperature of 70-80 K and density of 4x10^ cm~'^, embedded in surrounding lower 
density gas (n(H2) ~ 100 cm~^), and that the N159S clump shows nearly uniform-temperature 
as is consistent with no internal significant heat source. 

A number of attempts have been made to derive physical quantities in the N159 by using 
multiple CO emission lines and comparing them with line transfer calculations (Bolatto et al. 
2005; Nikolic et al. 2007; Pineda et al. 2008; Minamidani et al. 2008). Bolatto et al. (2005) 
used observations of the ^^CO J=4-3 and J=l-0 transitions at a 109" resolution and suggested 
that N159W has a cold component of 20 K and 10^ cm~^ and a hot component of 100 K and 
10^ cm~^. We present a comparison of the results, where we listed the most recent three works 
with high angular resolutions, in Table 4. The present work is distinguished from them in that 
it uses a combination of high excitation lines (the J=4-3 and J=3-2 lines) and optically thin 
^^CO lines at a high angular resolution of 45", nearly equal to that of Pineda et al. (2008). We 
also note that the signal to noise ratio of the present ^^CO J=4-3 data is high (minimum S/N 
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~ 8 at a resolution of 45") , helping to obtain better constraints on physical parameters. 

In N159W, Pineda et al. (2008) used ^^CO J=4-3, J=7-6, ^^CO J=4-3 and [Cl] ^P^- 
^Pq, ^P2-^Pi fine structure transitions observed by the NANTEN2 telescope and smoothed to a 
38" resolution, to derive a temperature of ~ 80 K and a density of ~ 10^ cm~^ toward N159W, 
with a PDR model. Minamidani et al. (2008) used ^^CO J=3-2, J=l-0 transitions and give 
constraints on temperature to be higher than 30 K and density in a large range from 3 x 10^ 
cm~'^ to 8 X 10^ cm~^. The present results show a fairly good agreement with the Pineda et al. 
(2008) 's results. We note that the Pineda et al. (2008) observed transitions up to ^^CO J=7-6 
line, higher than the present work, which may lead to sample denser gas than the present work. 
We also note that the present results have significantly improved the ranges of temperature 
and density as compared with Minamidani et al. (2008). A similar significant improvement is 
also found for N159E (Table 4). The present results for N159S are consistent with those of 
Minamidani et al. (2008) who used observations of ^^CO and ^^CO J=l-0 and ^^CO J=3-2 
transitions. This may refiect that N159S is a lower density cloud which can be well probed 
by the transitions lower than the J=3-2 line. A two-component model is proposed by Nikolic 
et al. (2007) which employed the ^^CO data up to the J=3 state and the CS data. We suggest 
that the model may be connected to the present idea of foreground self-absorption discussed 
above, while the physical parameters are yet poorly constrained. 

5. Comparisons with Star Formation Activities 

The N159 region is an active site of massive stare formation. The N159 itself is a group 
of Hll regions cataloged by Henize (1956). Bica et al. (1999); Bica et al. (2008) and Nakajima 
et al. (2005) presented catalogs of star clusters at optical and near infrared wavelengths. Radio 
studies revealed additional signs of young stars, including radio continuum thermal emission 
and H2O and OH maser emissions (Hunt & Whiteoak 1994; Caswell & Haynes 1981; Scalise & 
Braz 1981). More recently, Jones et al. (2005) made a study of young objects with the Spitzer 
IRAC, and presented a comprehensive view on young stars observable at near- to mid-infrared 
wavelengths. Figures 6(a) and 6(b) present distributions of these young objects, stars and 
clusters with nebulosity. Tables 5 and 6 list star clusters and Hll regions, respectively. The 
total luminosity of these objects are estimated to be ~ 7 x 10'^ Lq, most of which is emitted 
by the most massive O type member stars. In the following, we compare the molecular gas 
distribution with the young star clusters and the signs of star formation activities, including 
mid-infrared emission. 

5.1. Ha 

Figure 7 shows the distribution of the ^^CO J=4-3 integrated intensity superposed on 
an Ha image based on observations made with ESO Telescopes at the La Silla Observatory 
(programme ID 076.C-0888; processed and released by the ESO VOS/ADP group). The peak 
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of N159W has no Ha emission, but is associated with a bright, extended Ha emission centered 
at (RA, Dec) ~ (5^39""40^ -69° 46' 30") in the south. We also notice some faint Ha features 
around the molecular peak of N159W. The most significant Ha emission is extended around 
N159E with an extent of ~ 4' in RA and ~ 4' in Dec. This H II region has a dark lane extending 
in the east to west direction which coincides well with the N159E molecular clump. The peak 
molecular column densities toward N159W and N159E are estimated to be 2.2 x 10^^ cm~^ and 
1.4 X 10^^ cm~^, respectively. By using the relationship with visual extinctions (Imara & Blitz 
2007; Dobashi et al. 2008), we obtain extinctions of 20 - 110 mag and 18 - 72 mag for N159W 
and N159E, respectively, which are large enough to explain the dark lanes by the molecular 
gas. 

5.2. Mid-Infrared Data from Spitzer SAGE Program 

Mid-infrared measurements are one of the useful probes for investigating star formation 
activities, because the spectral energy distribution (SED) peak is located in this wavelength 
range (Calzetti et al. 2007). Spitzer SAGE program surveyed the entire LMC from 3.6 fim 
to 160 fim (Meixner et al. 2006), allowing a study of the star formation activity throughout 
the LMC (Whitney et al. 2008) and "molecular ridge" (Indebetouw et al. 2008). To identify 
embedded young stars in the observed region, we have inspected the Spitzer data at wavelengths 
shorter than 70 /im (Meixner et al. 2006), yeilding 5 mid- infrared peaks, seen as compact peaks 
at both 8 fim and 24 /xm. They are named from ^1 to #5 as listed in Table 7. Note that this 
naming is different from that used by Hunt & Whiteoak (1994). There is an embedded ultra- 
compact Hll region toward ^1 and another compact Hll region toward #5. Peaks #2/7^3 
and #4 are also associated with radio Hll regions. This indicates that all the five sources are 
associated with OB stars. 

Figure 8 shows a comparison of the ^^CO J=4-3 distribution with the Spitzer images at 
two wavelengths of 8/im and 24/im. The five Spitzer sources show fairly good spatial correlation 
with the 12 CO and are likely associated with the molecular gas. A more closer look shows that 
the molecular peak of N159W is significantly shifted from #1 by ~ 25" and from ^2 and #3 by 
~ 40" and that the molecular peak of N159E shows an extension nearly similar to the elongated 
feature associated with 7^5. #4 is on the other hand located between the two ^^CO clumps. 

Figure 9 shows a comparison between the Ha image and the Spitzer 24/im distribution. 
Most of the 24/im features seem well associated with the Ha. Two exceptions are in the dark 
lane toward N159E where Ha becomes weak and toward the peak of N159W. These reflect dust 
extinction as already noted (see Figure 7). We note the two small peaks of the 24/im emission 
at (RA, Dec) ~ (5^'39'^30^ -69° 47' 15") and (5M0™00^ -69° 47' 10"), small star forming regions 
in Figure 6 (a), show good correspondence. 

5. 3. Comparison of the CO Line Intensity Ratio with Star Formation 

Minamidani et al. (2008) showed good correlation between line ratio -R3_2/i-o and Ha 
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flux in clump-averaged values. Ha is however not a good tracer of star formation in high 
density regions because of the dust extinction. Figure 10 shows pixel by pixel comparisons 
of -R3_2/i-o which Minamidani et al. (2008) used 24yum flux. 24/im flux shows a fairly good 
correlation with a correlation coefficient of 0.62. Mid-infrared (MIR) emission is less sensitive 
to the dust extinction than shorter-wavelength emission, and is able to trace embedded young 
objects (Jones et al. 2005). We shall compare the CO line ratio with the MIR emission obtained 
by Spitzer SAGE observations. 

Figure 11 (a) shows the distribution of the integrated intensity for the line ratio, ^^CO 
J=4-3 to ^^CO J=l-0 (hereafter i?4_3/i_o) in the northern part of N159. Figure 11 (b) shows 
the same distribution, superposed on contours of the 24/im emission observed with Spitzer. 
In general, the ^^CO J=4-3 emission shows a good correlation with i?4_3/i_o. Figure 11 (b) 
reveals that MIR peak #5 and its associated extension at 24/im in N159E are in excellent 
agreement with the enhancement of i?4_3/i_o, suggesting that the heating by the 24/tm source 
leads to the enhanced i?4_3/i_o. 24/im peaks in N159W, #1 and #2/#3, are located toward 
the edge of the enhancement of -R4_3/i_o and a similar trend is found with _R3_2/i_o (Figure 
11(c)). The region of high line ratios with i?4_3/i_o greater than 1.0 is distributed in the ^^CO 
J=4-3 clump of N159W and shows no sign of star formation except for the three MIR sources. 
By considering that the high line ratios indicate higher density, this motivates us to suggest 
that the main part of the N159W cloud is prior to active star formation except for the three 
MIR sources located toward the eastern edge and a few faint Hll regions. Further detailed 
observations at higher resolutions of the N159W cloud should be of considerable importance in 
probing the initial conditions of massive star formation. 

Figure 11 (d) shows that OB stars and Herbig Ae/Be stars are mainly distributed be- 
tween the two molecular peaks, where the J=4-3 emission is weak. This suggests that the 
natal molecular gas of the young stars has been dissipated by the ionization and the dynamical 
effects of the stellar winds. 

We show the same comparison for N159S as shown in Figure 12. We find a slight 
enhancement of -R4_3/i_o but no significant hint of star formation. The marginal enhancement 
in i?4_3/i_o in the east is ascribed to a small Hll region (Minamidani et al. 2008). 

5.4- Star Formation in N159 and the Nature of the N159W Molecular Core 

N159 is one of the most active sites of star formation in the LMC (Kawamura et al. 
2009). The present study combined with the ^^CO J=3-2 study by Minamidani et al. (2008) 
and Pineda et al. (2008) provides a better insight into temperature and density in the N159 
molecular clumps by using the sub-mm transitions of CO at 10 pc resolution. 

The works on star clusters cited in Table 5 indicate that the N159 region is similar 
to the Galactic massive star forming regions, in terms of the level of activity and short time 
scale (Jones et al. 2005; Bica et al. 1999; Nakajima et al. 2005). The main properties of star 
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formation are summarized in Table 8, along with those with the rj Car region for comparison. 

We make a comparison of the N159 GMC and rj Car CMC, one of the most active 
massive star forming regions in the Galaxy. In Figure 13 the top-left diagram shows molecular 
distribution of N159W observed in the ^^CO J=3-2 emission line and the bottom-left shows 
that observed in the ^^CO J=l-0 emission line. The two distributions are converted into the 
same spatial resolution with the top diagram in the t] Car northern cloud (Yonekura et al. 
2005). While ^^CO J=l-0 and J=3-2 emission lines trace different excitation conditions, the 
distributions of these lines are similar with each other qualitatively as mentioned in Section 
3.1. We also show the intensity distributions along with the X and Y axes in top-right and 
bottom-right diagrams. 

Observational properties of star formation activity such as the number of O stars and 
IR luminosities and the properties of GMCs, such as sizes and masses, in the N159 region and 
r] Car region are summarized in Table 8. The total mass and size of the two cores are nearly 
the same, ~ IO^Mq and ~ 10 pc. The star clusters are mostly OB associations including 4-10 
03-05 stars over an area of ~ 100 pc^ and there is no sign of spatially confined super star 
cluster like R136. 

N159 is Type III in the GMC classification (Kawamura et al. 2009). The present work 
at higher resolution indicates that a GMC may be resolved into sub structures; clumps N159W, 
N159E and N159S. Such clumps may be classified in a similar way from Type I to III; N159W 
and N159E are Type III, and N159S Type I or II as has already been discussed by Minamidani 
et al. (2008). Future follow-up works at high resolutions will help to refine this classification. 

The present study has shown that the western part of the N159W core is not associated 
with any optical or infrared features, suggesting that it is at a stage prior to active star forma- 
tion. Further detailed observations at higher resolutions should be of considerable importance 
in probing the initial conditions of massive star formation. 

6. Summary 

We carried out ^^CO J=4-3 observations of the N159 region with NANTEN2 and ^^CO 
J=3-2 observations toward ^^CO J=3-2 peaks with ASTE. The main results are summarized 
as follows; 

1. The N159 GMC has been resolved into three prominent clumps, N159W, N159E, and 
N159S in the ^^CO J=4-3 emission line. N159W shows the highest intensities among the 
three in the CO transitions from J=l-0 to J=7-6. 

2. Molecular densities and temperatures have been derived toward the three peaks. Using a 
LVG analysis involving the ^^CO J=4-3 , ^'^CO J=3-2 emission lines in addition to the 
other CO lines published, we find that N159W and N159E have temperature of ~ 70 - 80 
K and density of ~ 3 x 10^ cm~^, and that N159S has temperature of ~ 30 K and density 
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of ~ 1.6 X 10^ cm~^. In order to explain lower line intensities than expected, we suggest 
that ^^CO J=l-0 and J =2-1 lines may be affected by self-absorption by foreground lower 
excitation gas. 

3. The ^^CO J=4-3 distribution is compared with Ha and mid- to far-infrared emission, a 
sign of embedded star formation, obtained with the Spitzer SAGE program. N159E shows 
a good coincidence with a dark lane of Ha and also with a 24 /im extended source. On 
the other hand, N159W is associated with three compact 24 /im sources and some small 
Ha features, although the ^^CO J=4-3 peak of N159W and its western part show no sign 
of star formation. 

4. A comparison between N159 and rj Car indicates that they show similar star formation 
activity and we do not see significant difference in physical parameters between these two 
massive star forming regions. 
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Table 1. Observed intensity peak properties 
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Table 2. Intensity and intensity ratio at clump peak 
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Table 3. The results of the LVG calculation in each clump 
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Note.- Calculated temperature and density from LVG simula- 
tion where value is minimum. iVco is calculated using av- 
eraged velocity dispersion derived by ^^CO J=3-2 observation of 
Minamidani et al. 2008. 



Table 4. Summary of excitation analysis 

N159W N159E N159S 

Paper ref. T^in n(Ri) Tkin n{H2) Tkin n{H2) 

(K) (10^ cm-3) (K) (10^ cm-3) (K) (10^ cm-^) 



This work 




72 tl 


4 0+0 


oq+26 
f5'3_20 


q 1 +1.9 

^■^ -0.6 


31 


+ 10 

-8 


1.6 t°o: 


Nikohc07 CDC 


(1) 


10 - 100 


> 1 






< 


20 


all 


HTC 


(2) 


>100 


>3 






> 


20 


<3 


MinamidaniOS 


(3) 


>30 


3 - 800 


>40 


1 - 300 


20 


- 60 


1 - 6 


PincdaOS 


(4) 


80 


10 - 100 













Reference — (1) Hot component of results from RADEX calculation by Nikolic et al. (2007). 
(2) Cold component of calculation by Nikolic et al. (2007). (3) LVG analysis by Minamidani 
et al. (2008). (4) Clumpy-PDR analysis by Pineda et al. (2008). 
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Table 5. Cluster catalog 



Positioii(J2000) Size Associated 



object 


R.A. 

(h m s) 


DEC. 

(0 , 


maj 
(') 


min 
(') 


ref 


Molecular clump 


LMC-N159H 


5 39 30 


-69 47 28 


0, 


.3 


0, 


.3 


1 


N159W 


LMC-N159K 


5 39 31 


-69 46 5 


0, 


.3 


0, 


.3 


1 


N159W 


LMC-N159J 


5 39 32 


-69 43 54 


0, 


.3 


0, 


.3 


1 


N159W 


LMC-N159A1 


5 39 39 


-69 45 52 


0, 


.3 


0, 


.3 


1 


N159W 


NGC2079 


5 39 39 


-69 46 22 


1 


.1 


0, 


.9 


1 


N159W 


NGC2078 


5 39 40 


-69 44 21 


1 


.2 


1 


.1 


1 


N159W 


LMC-N159F 


5 39 40 


-69 44 33 


0, 


.8 


0, 


.7 


1 


N159W 


LMC-N159A2 


5 39 40 


-69 46 56 


0, 


.3 


0, 


.3 


1 


N159W 


NGC2084W 


5 39 53 


-69 45 42 


1 


.5 


1 


.3 


1 


N159E 


BSDL2727 


5 39 54 


-69 43 22 


0, 


.4 


0, 


.4 


1 


- 


NGC2083 


5 39 59 


-69 44 9 


2, 


.0 


1 


.8 


1 


N159E 


LMC-N159E 


5 40 2 


-69 47 12 


0, 


.6 


0, 


.6 


1 


- 


NGC2084e 


5 40 6 


-69 45 51 


1, 


.7 


1, 


,2 


1 


N159E 


BSDL2755 


5 40 19 


-69 46 26 


0, 


.8 





.6 


1 




N159-Y4 


5 40 4 


-69 45 8 


1, 


.1 


1, 


.1 


2 


N159E 


LMC-N159L 


5 40 2 


-69 49 10 


1, 


.4 


1, 


.1 


1 


N159S 


LMC-N173 


5 40 20 


-69 53 8 


0, 


.9 


0, 


.9 


1 


N159S 


LMC-DEM279 


5 40 26 


-69 50 22 


1, 


.2 


1, 


.2 


1 




N159S-Y1 


5 40 13 


-69 49 57 


0, 


.6 


0, 


.6 


2 


N159S 



Reference — (1) Bica et al. 1999; (2) Nakajima et al. 2005. 
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Table 6. HII regions 





Position(J2000) 


Size 




Ionizing 




Associated 


object 


R.A. 


DEC. 


maj 


min 


ref 


Source 


ref 


Molecular clump 




(h m s) 


(° ' ") 


(") 


(") 










LML-N159H 


5 39 30 


-69 47 28 






1 






AT "I P" ATS 7 

N159W 


T T\ /T/^ Ml CAT/' 


5 39 31 


-69 46 5 


25 


22 


1 






ATI c mn7 

W159W 


LMC-N159J 


5 39 32 


-69 43 54 






1 






AT 1 P" ATX 7 

N159W 


T ~\ If Ml cn A 1 

LMC-iNl59Al 


5 39 39 


-69 45 52 






/ 






ATI C ATlt 7 

M159W 


T 1\ /T/^ M1 r/^ A 

LML-N159A 


5 39 39 


-69 46 22 


56 


62 


1 


05V+07V 


3 


AT 1 X" AT T 7 

JN159W 


T 1\ /T/^ M1 r/^TT^ 

LML-N159i' 


5 39 40 


-69 44 33 


48 


47 


1 






AT "I P" ATt 7 

N159W 


T T\ If Ml A o 


5 39 40 


-69 46 56 






I 


(J9V 


3 


ATI C ATI 7 

M159W 


T ~\ If ATI f i'\ f ^ 


5 39 53 


-69 45 42 


140 


107 


1 






N159E 


LMC-M159D 


5 39 59 


-69 44 9 


131 


97 


1 






AT-1 CATT' 

N159E 


T T\ If ATI P'ATT' 


5 40 2 


-69 47 12 


35 


35 


1 








T H If AT "I TAT 

LMC-N159L 


5 40 2 


-69 49 10 


62 


28 


1 






ATi rAO 

N159h 


LMC-N173 


5 40 20 


-69 53 8 


49 


54 


1 


- 


- 


N159S 


LMC-DEM279 


5 40 26 


-69 50 22 


65 


65 


4 








N159AN 


5 39 37 


-69 45 26 


3.6 


3.1 


5 


2x04V/05V 


6,7 


N159W 


Papillon 


5 40 5 


-69 44 38 


3.1 


2.7 


8 


04V/05V 


6 


N159E 


Reference — (1) 


Henize 


1956; (2) 


Heydari- 


-Malay( 


2ri & 


Tester 1982; 


(3) 





Deharveng & Caplan 1992; (4) Davies et al. 1976; (5) Hunt & Whiteoak 1994; 
(6) Marti'n-Hernandez et al. 2005; (7) Indebetouw et al. 2003; (8) Heydari- 
Malayeri et al. 1999. 



Table 7. Mid-infrared peaks 





Peak position (J2000) 


Associated 


Associated 


IR 


Peak num. 


R.A. 


DEC. 


Molecular clump 


H II region 


Luminosity 




(h m s) 


(0 , 






(xIO^Lq) 


#1 


5 39 38 


-69 45 26 


N159W 


N159AN, HWt#4 


3.7 


#2 


5 39 37 


-69 46 8 


N159W 


N159A, HW#5 


1.7 


#3 


5 39 42 


-69 46 11 


N159W 


N159A, HW#5 


1.3 


#4 


5 39 52 


-69 45 17 




HW#2/#3 


1.2 


#5 


5 40 05 


-69 44 38 


N159E 


Papillon nebula, IIW#1 


3.4 



t HW means Hunt & Whiteoak 1994. 
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Table 8. Comparisons with r] Carinac GMC 





N159 region 


rei 


r\ Car region 


ref 


Massive star 


3x 04/05 


1 


5x03+7^Car 


2 


iNum. 01 U star 


> 55 


3 


70 


4 


oize 01 LjIVLL^ 










(pcxpc) 


ioUxyu 


5 


i5UX iUU 


b 


Mass of GMC 










(xIO^Mq) 


5.3 


5 


3.4 


6 


Total IR luminosity 










(xlO^Lo) 


~0.9 


7 


1.2 


8 



Reference — (1) Estimated from N^yc derived by Indebetouw et al. 
2003 and Martin-Hernandez et al. 2005; (2) Smith 2006; (3) Nakajima 
et al. (2005); (4) Simbad astronomical database; (5) A/vir from ^^CO 
J=3-2 observation by Minamidani et al. 2008; (6) Mco from ^^CO 
J=l-0 observation by Yonekura et al. 2005; (7) Jones et al. 2005; (8) 
Smith & Brooks 2007. 
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12CO/=1-0 l2co/=3-2 12C0/=4-3 




I I I I I I I I I 1 I I I I I I I I I I I I I I I I I I L 

5h 40m 30s 40m 00s 39m 30s 5h 40m 30s 40m 00s 39m 30s 5h 40m 30s 40m 00s 39m 30s 



Right Ascension (J2000) 

Fig. 1. Integrated intensity map of ^^CO J=l-0 (left), ^^co J=3-2 (middle) and ^^CO J=4-3 
(right). Lowest contour levels are 3cr (11, 5, 5 K km s^^for ^^co J^l-Q, ^'^CO J=3-2 
and ^^CO J=4-3 respectively) level of each observation and increments are 10 K km s^^. 
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(a) N159W 



4.0 



I- 0.0 



. (da, d5)=(20", 20"). 


(0", 20"). 


(-20", 20"). 


(20", 0") 


nJ] (0", 0"). 


(-20", 0") 


(20", -20") 


(0", -20") 


(-20", -20") 



220 240 260 

Vlsr (km/s) 
(b)N159E 



2.0 



I- 0.0 



. (da, d8)=(20", 20"). 


^^^^ 


(0", 20"). 


( 


-20", 20"). 




(20", 0") 




\ (0", 0"). 




(-20", 0"). 


(20", -20") 




(0", -20") 


(-20", -20") 



220 240 260 

Vlsr (km/s) 

Fig. 2. Profile map of ^^CO J=3-2. (a) N159W (center position is 5'^39"^36.8% -69° 45' 32" at 
J2000 coordinate), (b) N159E (5'^40™8.7% -69° 44' 34"), and (c) N159S (5h40'"5% -69° 50' 34"). 
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(c)N159S 



(da, d5)=(20", 20") 


(0", 20") 

■n jJ\ 


(-20", 20") 

■ « ■ 


(20", 0") 


(0", 0") 

. A ^ 


(-20", 0") 


(20", -20") 


(0", -20") 


(-20", -20") 



220 240 260 

Vlsr (km/s) 

Fig. 2. continued 
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220 230 240 250 220 230 240 250 220 230 240 250 



VIsr [km/s] 

Fig. 3. i^co ,7=1-0, i^co J=3-2, ^^CO J=4-3 and ^^CO J=3-2 line spectra 
from peak position of N159W (left), N159E (middle) and N159S (right). All spec- 
tra were convolved to 45" resolution. ^"^CO J=3-2 spectra were multiplied by 5. 
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One component 



High excitation component 



100 



c 




100 



10 























- (c) N159E 

"i 1 1 1 1 1 


^ (d)N159E - 



100 



c 



10 





3 4 5 

log n(H2) [cm"^] 



3 4 5 

log n(H2) [cm"^] 



Fig. 4. Observed line intensity ratio overlaid with distributions of the reduced value. Green 
lines shows a subset of line ratios. Alphabets correspond column (5) of Table 2; (A) ^^CO 
J=4~3/i2CO ,7=1-0, (B) i^co J=4-3/i2CO ,7=2-1, (C) ^^co 7=4-3/i2CO 7=3-2, (D) ^^co 
7=7-6/i2CO J=4-3, (E) ^^co 7=4-3/i3CO J=l-0, (F) ^^co J=4-3/i3CO J=2-l, (G) ^^co 
J=4-3/i3CO J=3-2, (H) i^co 7=4-3/1^00 J=4-3. Red contour means threshold value of 
X^. Left diagrams are the results of LVG simulation using all observed lines. Right diagrams 
are the result using ^^CO J=3-2 and higher excited rotational transitions and all ^'^00 lines. 
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CO 
CO 


c 



00 



CO 
CO 


c 



00 



8 

6 

4 

2 


10 

8 

6 

4 

2 



N159W 




-I 1 1- 



T 1 r- 



N159E 




-I 1 1- 



N159S 




1-0 2-1 3-2 4-3 7-6 

Transition of CO 

Fig. 5. Plots show calculated intensity from LVG with best fit parameters of high ex- 
citation component (green line) and observed line intensities (mark with error bar). 
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^ I I I I I I I I I ~'< I I I I I I I I I I I I I I 

5h40m15s 40m00s 39m45s 39m30s 39m15s 5h40m15s 40m00s 39m45s 39m30s 39m15s 

R.A. (J2000) R.A. (J2000) 



Fig. 6. (a) Distributions of massive stars and young stellar objects (*1 Feast et al. 1960, 
*2 Walborn 1970, *3 Gatley et al. 1981, *4 Deliarveng & Caplan 1992, *5 Cowley et al. 
1995, *6 Lazendic et al. 2002, *7 Jones et al. 2005, *8 Martin-Hernandez et al. 2005, 
*9 Nakajima et al. 2005, *10 Tester et al. 2007, *11 Tester et al. 2006 overlaid with 
Ha image. (b) Distributions of cluster with nebulosity cataloged by Bica et al. (1999). 
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5h40m00s 39m30s 

R.A. (J2000) 



Fig. 7. ^^CO J=4-3 integrated intensity map overlaid with Ha image. 
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Fig. 8. Contour of ^^CO 7=4-3 integrated intensity overlaid with 8/im (a) and 24/im 
(b) image from Spitzer/SAGE. Properties of marked peak are listed in Table 7. 
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Fig. 9. Diagram shows Ha image and contour of 24/im flux. Contour levels are 200, 400, 800, 1600 
MJy/sr. 
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2.5 
2.0 
O 1.5 

or 

1.0 
0.5 
0.0 



c.c.= 0.62 



10 100 

24 i^m flux [MJy] 



Fig. 10. Pixel by pixel 
with 24 /im flux, 



comparison 
where 



(20" 
the 



grid) of 
correlation 



line ratio 
coefficient 
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_] I I '"^ IT I I I I I 

5h40m00s 39m30s 5h40m00s 39m30s 

R.A. (J2000) R.A. (J2000) 



Fig. 11. (a) Image of integrated intensity ratio R4-3/1-0 and contour of ^^CO J=4-3 integrated in- 
tensity in northern part. (b) Contour of 24 /im flux superposed on R4„3/i„o and marked num- 
ber show 24/im intensity peak position, (c) R3_2/i_o and contour of 24 ^m. flux, (d) Positions 
of OB stars and Herbig Ae/Be stars derived by Nakajima et al. (2005) overlaid with R4_3/i_o. 
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O OB stars + Herbig Ae/Be stars 




5h40m20s 40m00s 5h40m20s 40m00s 



R.A. (J2000) R.A. (J2000) 

Fig. 12. Same as Figure 11 but in southern part. 
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Distance from peak (pc) 



Fig. 13. Molecular gas distributions of N159 (top) and rj Car, northern cloud (Yonekura et al. 
2005). The map of rj Carinae region was smoothed to same spacial resolution of N159 
(5pc). Contour levels are 10%, 20%, 90% of peak integrated intensity. Intensity dis- 

tribution in the directions of X, Y slit in the left figure are displayed in right figure. 



36 



